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It  is  shown  how  to  single  out  a  particular  diagonal  of  a  biinfinite  band 
matrix  A  as  its  main  diagonal,  using  the  decomposition  of  the  solution  set 
of  Ax^  =  (1  into  those  which  are  bounded  at  00  and  those  which  are  bounded 
at  -*>.  As  an  application,  it  is  proved  that  the  inverse  of  the  coefficient 
matrix  for  the  system  satisfied  by  the  B-splines  coefficients  of  the  cubic 
spline  interpolant  at  knots  is  checkerboard  and  that,  under  certain  assump¬ 
tions,  the  local  mesh  ratio  must  be  bounded. 


W 


AMS  (MOS)  Subject  Classifications:  Primary  47B37,  41A15,  47B39. 

Secondary  15A09,  39A12,  39A70,  47B37. 

Key  Words:  Biinfinite,  Banded,  Matrix,  Exponential  Decay,  Dichotomy,  Spline 
Interpolation,  Cubic,  Mesh  Ratio  Restriction. 

Work  Unit  Number  3  -  Numerical  Analysis 


/ 


Sponsored  in  part  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024 . 
Some  of  the  work  was  done  while  the  author  was  a  very  contented  and  grateful 
guest  of  the  SFB  72  at  the  University  of  Bonn,  Germany  in  the  summer  of  1979.  ^  [H 


SIGNIFICANCE  AND  EXPLANATION 


Spline  approximation  is  often  most  effective  when  the  breakpoint  (knot) 
sequence  can  be  chosen  suitably  nonuniform.  At  the  same  time,  standard  spline 
approximation  schemes  (such  as  least-squares  approximation  by  splines)  are  so 
far  only  known  to  be  bounded  as  long  as  the  breakpoint  sequence  is  almost  uni¬ 
form.  Any  such  bound  is  obtained  (explicitly  or  implicitly)  in  terms  of  a 
bound  on  the  inverse  of  certain  matrices  which  are  banded.  Any  attempt  at 
establishing  bounds  for  more  general  breakpoint  sequences  must  therefore  come 
to  grips  with  the  inverses  of  these  band  matrices.  The  hope  is  the  Demko's 
discovery  of  the  exponential  decay  of  band  matrix  inverses  will  lead  event¬ 
ually  to  those  desired  bounds. 

In  the  present  report,  a  specific  question  concerning  the  boundedness 
of  cubic  spline  interpolation  at  breakpoints  in  terms  of  the  local  mesh  ratio 
leads  to  a  description  of  the  inverse  of  a  biinfinite  band  matrix  A  in  terms 
of  the  behavior  of  the  solutions  of  the  homogeneous  problem  Ax  =  0.  This 
description  should  be  of  use  in  the  analysis  of  other  spline  approximation 
schemes  as  well. 


»T.tFSSI0N  (or  y 

NTIS 

DDC 

UNM1N0UW 

JUSUFICA7 

White  Section 

Butt  Section  □ 
DED  □ 

BY 

DISTfflBUT! 

Dist.  A\ 

ON/WAjUWUn  EBB 

The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


WHAT  IS  THE  MAIN  DIAGONAL  OF  A  BI INFINITE  BAND  MATRIX? 


Carl  de  Boor 

* 

INTRODUCTION 

'  The  study  of  approximation  by  splines  on  a  biinfinite  knot 

sequence  leads  to  linear  systems 
Ax  =  b 

with  a  banded  biinfinite  coefficient  matrix  A.  Questions  as  to 
the  existence  of  A  \  its  boundedness  or  its  possible  checker¬ 
board  nature  need  to  be  answered,  and  these,  in  turn,  raise  the 
question  of  which  diagonal  of  A  may  be  the  main  diagonal. 

For  example,  it  is  well  known  that  the  inverse  of  a  finite 
totally  positive  matrix  A  is  checkerboard,  i.e., 

(-)1+^A  1(i,j)  >_0,  all  i, j 

•  -1 

and,  in  particular,  the  entries  of  the  main  diagonal  of  A  are 

positive.  One  would  expect  the  same  statement  to  be  true  when 

A  is  biinfinite,  but  it  is  not  clear  a  priori  which  diagonals 

of  A  ^  will  be  positive  and  which  negative. 

Again,  in  one  approach  toward  proving  that  the  inverse  of 

a  biinfinite  totally  positive  matrix  A  is  checkerboard,  one 

would  try  to  show  that  the  inverse  is  approximated  in  some 

pointwise  sense  by  the  inverse  of  finite  sections  of  A,  whose 

checkerboard  pattern  is  then  known.  Now,  as  we  will  show 

? 

i 
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below,  if  there  is  any  convergence  at  all,  then  these  finite 
sections  are  necessarily  principal  with  respect  to  one  fixed 
diagonal  of  A  .  That  diagonal  then  has  earned  the  epithet 
"main" . 

Here  is  an  outline.  In  Section  1,  even-order  spline  in¬ 
terpolation  at  knots  is  discussed,  since  I  was  led  to  wonder 
about  the  main  diagonals  of  biinfinite  matrices  because  of  an 
argument  with  C.A.  Micchelli  concerning  mesh  ratio  restrict¬ 
ions  for  that  scheme.  Section  2  contains  a  discussion  of  bi¬ 
infinite  band  matrices,  in  particular  some  propositions  re¬ 
garding  existence  and  character  of  their  inverses.  In  the  last 
section,  some  of  these  results  are  applied  to  cubic  spline  in¬ 
terpolation  at  a  biinfinite  knot  sequence,  giving  me  an  oppor¬ 
tunity  to  correct  a  mistake  in  [4] . 

1.  MESH  RATIO  RESTRICTIONS 
IN  EVEN-ORDER  SPLINE  INTERPOLATION  AT  KNOTS 

Let  m$k  t  be  the  normed  linear  space  of  bounded  splines 
of  order  k  with  knot  sequence  t  =  (t^)  and  the  sup-norm. 

We  take  the  knot  sequence  t  to  be  biinfinite  and  strictly 
increasing.  Also,  let  t  be  a  strictly  increasing  biinfinite 
sequence.  The  problem  is  to  determine  for  given  a  e  an  f 
in  m$.  with  f(x.)  *  a.  ,  all  i.  I  call  this  interpolation 

K  /  u  1  1 

problem  correct  (others  have  called  it  "poised") if  it  has  ex¬ 
actly  one  solution  for  every  afci,,,  i.e.,  if 

R  }  mSk,t-^  *«!  f|x  s=  (f<Ti>>-” 

is  invertible.  Note  that  R*1  is  necessarily  bounded  if  it 
exists. 

This  interpolation  problem  has  received  particular  atten¬ 
tion  in  the  special  case  of  even-order  spline  interpolation  at 
knots,  i.e., 

x_  »  t  and  k  even  . 

See,  e.g.,  [6]  for  a  recent  survey.  I  showed  in  my  talk  1 4 J  at 
the  last  Bonn  conference  that  this  interpolation  problem  is 
correct  in  case  the  global  mesh  ratio 


is  finite.  I  also  stated  there  without  proof  that,  for  a  cor¬ 
rect  problem,  the  local  mesh  ratio 


m.  :=  sup  At. /At. 

would  have  to  be  finite,  since  it  is  possible  to  bound  this 
mesh  ratio  in  terms  of  ||  R  .  It  was  this  claim  which 
Micchelli  doubted  when  we  discussed  various  possible  suffi¬ 
cient  conditions  for  the  correctness  of  the  interpolation 
problem  last  summer.  Now,  my  claim  was  based  on  the  corre¬ 
sponding  result  in  case  of  a  finite  knot  sequence,  in  [3], 

Here  is  an  adaptation  of  the  argument  there  to  the  present  bi¬ 
infinite  context. 

Supposing  the  problem  correct,  write  the  interpolant  R 
to  the  particular  data  sequence  £  :=  ((-)1)  in  terms  of  the 
normalized  B-splines  of  order  k  for  t  , 

R-1f  =  £  a .N .  , 

j  3  3,k 

Then,  from  [2], 


I  R_1iH  oo  i 


|a| 


00 


for  some  positive  constant  independent  of  t  . 

any  particular  i  , 

-1  -1  ai-oti 

=  tvw*  £  -  <*  £>  ’ (5)  -  ztT^ 

Since  (N .  .  . )  forms  a  partition  of  unity  and  N . 

j  »*”■*•  3 

its  support  in  (t^t^^)  while  Ult.,t.+1)  , 

that 


Now,  for 


,k-l 


U> 


,k-l  has 
this  implies 


min 

i-k+2<j<i 


j-t-k-1 


-t . 


TtT- 


i  “k 


dJR_1S 


and  it  was  from  this  inequality  that  I  had  drawn  a  bound  for 
~1|j  .  But,  actually  this  inequality  is  not 
strong  enough  for  such  a  conclusion  in  case,  e.g.,  tj+k_^-t^ 
is  a  decreasing  function  of  j  .  The  desired  conclusion  can  be 


m.  in  terms  of  Jr 
t 


drawn,  though,  if  we  are  certain  that  the 
alternate  in  the  right  way. 


‘j 


alternate,  and 


PROPOSITION  1.  If  the  B-spline  coefficients  a  of  the  inter¬ 
polant  R_1f  to  the  data  sequence  f  ■  ((-)i)  alternate  in 
such  a  wav  that  >  0.  all  i.  then  mfc  <_  DjJ|R  II  • 


Proof.  Under  this  additional  assumption,  we  have 

(-) £  0  for  j  =  i  and  also  for  j  =  i-k+2  (since 

k  is  even) ,  hence  from  the  above 


i-1 

I 


j=i-k+3 


(-) 


1*1  Wi 

tj+k-rtj 


,k-i 


<C> 


t  2ll“L/min(tj+k-i_tj  :  i-k+3  1  3  <  i-1) 
and  therefore 


Ati-1  +  Ati+1 
Ati 


<_  min 
i-k+3< j<i-l 


ti-t-k-rti 


Ati 
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which  does  imply  the  desired  result. 


This  leaves  open  the  question  as  to  when  we  can  expect 
the  alternation  assumption  to  hold,  in  a  finite-dimensional 
situation,  the  alternation  is  immediate  because  of  the  total 
positivity  of  the  coefficient  matrix  of  the  linear  system 


(-)1 


all  i 


But  the  sense  of  the  alternation  depends  on  just  what  the  range 
of  i  and  j  here  is.  Specifically,  if  i  *  I+l,...,I+n, 
j  ■  J+l,  . ..,  J+n,  then  (assuming  that  J  <  I  <  J+k) 


(-)I+j  aJ+j  >  0  ,  all  j  . 

In  the  biinfinite  case,  it  is  not  even  clear  that  a  al¬ 
ternates.  In  order  to  investigate  this  question  further,  we 
now  turn  to  an  analysis  of  biinfinite  band  matrices. 


2.  THE  INVERSE  OF  A  BI INFINITE  BAND  MATRIX 


in  order  to  give  our  definition  of  main  diagonal,  we  need 
notation  to  describe  sections  of  biinfinite  matrices  and  their 
relationship  to  the  linear  map  represented  by  that  matrix. 

Let  I ,  J  be  integer  intervals.  Then 

AI,J  *"  A|lxj  “  >iti#  jcj 

denotes  the  corresponding  section  of  the  biinfinite  matrix  A. 
He  can  think  of  Aj  j  simply  as  a  |l|x|j|  matrix.  But,  A 1J 
also  describes  the  nontrivial  part  of  the  linear  map 
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f 


* 


i- 


with 


(Pja) (i) 


i  e  i 

i  ig¬ 


nore  precisely,  AT  is  the  matrix  representation 

i.  ,  J 

spect  to  the  canonical  basis)  of  the  linear  map 


(with  re- 


Pj  (A 


ran  P 


) 


and  we  will  not  distinguish  between  these  two. 

Here  and  below,  we  will  use  the  alternative  notation  a(i) 
for  the  i-th  entry  of  the  sequence  a  ,  which  is  consonant 
with  the  notation  A(i,j)  for  the  (i,j)-entry  of  the  matrix  A. 


DEFINITION.  The  bounded  and  boundedly  invertible  biinfinite 
matrix  A  (as  a  map  on  i^fE)  «  say)  has  its  r-th  diagonal  as 
main  diagonal  ;=  the  matrices  (AT  _  )  ^  converge  finitely 

to  A  as  I XL  ,  i.e.,  A£  I+r  is  invertible  for  all  suf¬ 
ficiently  large  finite  intervals  I  and 

A_1(i,j)  =  lim  (A  T  )_1(i,j)  ,  all  i,j  . 

I-*-  7L  1»'L+r 

Here,  we  have  used  the  abbreviation 
I+r  :=  {i+r  :  iei}  . 

For  example,  suppose  that  both  A  and  A-1  are  upper 
triangular  in  the  usual  meaning  of  the  word,  i.e.  , 

A(i,j)  =  A_1(i,j)  =  0  for  i  >  j  . 

Then 

j  -1 

6..  -  Z  A(i,k) A  1 (k, j ) 

13  k=i 

showing  that  (A  1)I  1  ■  (Aj  j)-*  .  In  this  case  then,  diago¬ 
nal  0  is  the  main  diagonal  of  A  {as  usual!).  But  now  let 
E  be  the  map  or  matrix  of  the  left  shift. 


(Ea) (i)  :=  a(i+l)  ,  all  i  , 

and  let  re  Z  .  Then  ECA  is  also  invertible,  with  inverse 
-1  -r 

A  E  .  But  now 

(ErA)  (i,  j)  *>  A(i+r,j)  ,  (ErA)'1(i,j)  «A_1{i,j+r) 


hence  now 


( (ErA) 


i1"1  ■  <Ai,irl 


'  <ErA1'l,I-r 


"  i-r , I  I , I  I, I  '  '  I , I-r 

In  other  words,  ErA  has  diagonal  r  as  its  main  diagonal 

j-  -  1 

(while  (E  A)  has  diagonal  -r  as  its  main  diagonal) . 

For  matrices  which  are  not  triangular  with  triangular  in¬ 
verse,  it  is  much  more  difficult  to  ascertain  whether  or  not 
they  even  have  a  main  diagonal,  let  alone  which  diagonal  it 
might  be.  We  now  discuss  this  question  in  the  context  of  band¬ 
ed  matrices. 


DEFINITION.  We  say  that  the  biinfinite  matrix  A  is  m-banded 
if 

(1)  A (i ,  j  )  =  0  for  j  j.  [i,i+m] 

(2)  A(i,i)A(i,i+m)  ^  0  ,  all  i  . 

Thus  an  m-banded  matrix  has  at  most  m+1  nonzero  bands. 

If  A  is  m-banded,  then,  for  any  r  ,  ErA  also  has  just  m+1 
possibly  nonzero  contiguous  bands  and  so  could,  with  justifi¬ 
cation,  also  be  called  m-banded.  But  we  will  use  the  term  "m- 
banded"  only  as  described  in  order  to  suppress  an  additional 
inessential  parameter.  The  other  assumption,  viz.  the  nonva¬ 
nishing  of  the  first  and  last  band,  is  a  nontrivial  one.  It 
makes  certain  statements  simpler  and  is  satisfied  in  the  case 
of  spline  interpolation  at  knots. 

An  m-banded  matrix  A  gives  rise  to  a  linear  map  on  S? 
which  we  will  identify  with  A  .  This  map  has  an  m-dimensional 
null space  or  kernel, 

tl:=  \  .  '  (felP  :  Af  =  0} 

In  particular, 

(3)  for  every  i,  — ►Kln  :  ft — *£|  [i+1  |i+mj  is  1-1  and  onto 

because  of  (2).  Conversely,  if  Xl  is  an  m-dimensional  subspace 
of  I?  for  which  (3)  holds,  then  there  is,  up  to  left  multi¬ 
plication  by  a  diagonal  matrix,  exactly  one  m-banded  matrix 
having  tl  as  its  nullspace. 

We  also  introduce  two  subspaces  of  1(1  , 

<l+  :*  {f  cIR.  :  Tim  f  (i)  <  »>}  ,  m+  :=  dim  #+ 

i+oo  ^ 

1 c  !-  {f«H  s  Ilm  f(i)  <  »}  ,  m“  :=*=  dim  ft 


0 


« 
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From  now  on,  we  assume  that  A  is  m-banded  and  bounded 
on  Eb  ,  Then  the  i-th  column  A  1  ( ■  , i )  of  its  inverse,  if  it 
exists,  would  solve  the  linear  system 

Ax  -  («i3>  • 

The  following  proposition  is  therefore  a  first  step  toward 
understanding  the  inverse  of  an  m-banded  matrix. 


PROPOSITION  2.  For  all  i,  there  exists  exactly  one  L.e£ 

-  - - — _ —  -j. 

such  that  AL .  =  (3.  .)  if  and  only  if  tl  =  >1  ©ft.  . 
-  _v  _  ^  - i - 

Proof .  Since  tl  rtf l  is  the  kernel  of  AL  ,  there  is  at 

'  cc 

most  one  solution  (for  any  particular  l)  if  and  only  if 

tfnn+  =  {0}  .  Hence  it  is  sufficient  to  prove  that,  given 
uniqueness,  there  is  a  solution  for  every  i  iff  m  +m+  =  m  . 


For  this,  note  that  L.  €.  1 

—l  = 


satisfies  AL.  =  (S..)  iff 


(4)  L^j)  = 


L7 < j )  ,  j<i+m 
L^(j>  ,  j>i 


,  with 


l7  £  tl 
Lte  H+ 


i+m 

(5)  1  A  ( i  ,  j  )  L .  ( j  )  =  1  . 

j  =  i 

*•*  + 

Here,  L.,  L.  is  the  extension  of 


-i  j  [i,i-t-m-l] 


— i| [i+1  i+m]  '  respectively,  to  an  element  of  tl . 

Now  to  see  that  m~+m+  =  m  implies  existence  of  L^  , 

note  that  (4)  and  (5)  constitute  a  linear  system 

O' 

<6>  VLlu.iW  '  ; 


in  the  m+1  unknowns  L^(i+m)  ,  with  the  first 

m-m-  homogeneous  conditions  ensuring  that  the  extension 

of  L. j  , .  . .  to  an  element  of  Yl  lies  in  ft”,  and  the 

—l  j  [i+i+m-l J 

next  m-m  homogeneous  conditions  ensuring  that  the  extension 

I,t  of  1.1,.,.  ,  to  an  element  of  H  lies  in  ft+,  and  the 

—i  —i I li+i , i+mj 

last,  the  only  inhomogeneous,  condition  being  (5).  But  if  now 
m~+m+  =  m  ,  then  (6)  has  as  many  equations  as  unknowns  and,  as 
we  already  know  that  it  has  at  most  one  solution,  the  existence 
of  a  solution  follows. 

Conversely,  assuming  the  existence  of  a  solution  for  every 
i  ,  consider  the  maps 
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♦*  :  rnm+1- 


>u 


m  * 

I  a  .  L  .  • 

3-1+3 


j  =  0 


with  *  standing  for  +  or 


,  and  m 


+ 

m  +m 


Then 


* 

dim  ker  $ 


m+1  -  dim  ran 


* 


_  * 
>  m+l-m 


*  =  + 


Consequently,  there  exists 


a  eker 


+ 

$ 


A  ker 


4>  \{0}. 


M 


m 

l 

3=0 


a  .  L .  ,  . 
D-i+3 


¥  0 


For  this  a, 


since  (L^)  is  obviously  linearly  independent.  On  the  other 
hand,  since 


we  find 


— i+j  ^  '  ^or  s  <  ^+m 
— i+j  ^  '  for  s  > 


j-0 , . . . ,m  , 


M  .  fX°  aj— i+j  °n  h“’i+mt 
U;  ajiii+2  on 

and  therefore,  by  choice  of  a  ,  M(s)  =  0  for  s  <  i+m  and 
s  >  i+m  .  This  implies  in  >_  m  ,  and  therefore,  since  by  as¬ 
sumption  >tlf  =  {0},  i.e.,  m++m  m,  the  conclusion 
m  +  m  *  m  follows. 


Next,  we  characterize  bounded  invertibility  of  a  bounded 
m-banded  matrix  A  in  terms  of  Yl . 

PROPOSITION  3.  A  jis  boundedly  invertible  if  and  only  if 

£i)  ^-=  ^  ©^  ,  (ii)  the  elements  of  tL  and  fl+  decay  expo- 

nentially  and,  (iii)  for  each  i,  the  matrix  in  (6)  can  be 

chosen  so  that  sup  |b-1|  <  »  . 

i  1 

Concerning  the  exponential  decay,  I  had  proved  at  the  last 
Bonn  conference  that,  for  a  bounded  and  boundedly  invertible 
m-banded  matrix  A  on  some  with  p  <  »,  there  exists 

const  so  that,  for  all  ftiC  and  all  i  , 

ll£(r)llp  1  const  Ar  I|  f (0)||P  ,  r-1,2,3,... 

with 

A  -Irr^l  '  K  8“  liA|1||A’111 


» 


9 


I 


i 


-0 


and 


f(r)  •=  fl 

—  '  — | [i+rn, i+ (r+l)m-l]  ' 

and  an  analogous  statement  for  n+  .  The  exponential  decay 
mentioned  in  condition  (ii)  of  the  proposition  is  meant  in 
this  sense. 

The  proof  of  the  necessity  of  the  exponential  decay  rests 
on  Demko's  [5]  nice  idea.  As  a  footnote,  I  would  like  to  re¬ 
cord  here  that,  in  response  to  my  talk,  S.Demko,  at  the  end  of 
his  talk  at  the  present  conference,  made  the  point  that  he  had 
been  materially  helped  by  a  referee's  report  authored,  as  it 
happens,  by  T.  Lucas. 

Finally,  the  proof  of  the  necessity  of  the  last  condition 
is  a  bit  tricky. 

On  the  other  hand,  the  sufficiency  of  the  three  conditions 
is  immediate  since  they  insure  that  the  L^ 1 s  ,  constructed  in 
Proposition  2  on  the  strength  of  (i) ,  are  t^-bounded  uniformly 
in  i  and  decay  exponentially,  hence  the  matrix  [...,L^,...] 

is  a  bounded  map  on  l  for  every  p  ,  etc. 

P 

The  complete  proof  of  Proposition  3  (and  of  Proposition  4 
to  follow)  can  be  found  in  [7] .  Finally,  we  state  a  necessary 
and  sufficient  condition  for  such  a  bounded  m-banded  matrix  A 
to  have  a  main  diagonal. 


PROPOSITION  4.  A  has  diagonal  r  as  its  main  diagonal  if  and 

+ 


only  if  (i)  A  is  boundedly  invertible,  (ii)  r  =  m  _ 

exists  a  positive  const  so  that,  for  all  large  n 


(iii)  there 


for  all  f  c  fl. 


and 


| l-n,-n+m+-l]"  — 


>  const 


-  [-n,-n+m] 1 


for_all  fett"  ,  111)  [n+m++1,n+m]«  i  const  li|[n,n+m]' 


In  effect,  under  these  assumptions,  we  can  construct  the 
{I)  of  (Ax  I+r)_1  for  all  large  I 
tion  of  the  corresponding  column  of  A-'*'  , 


column  L.' 

— l 


as  a  modifica- 


(I) 


"  h  ~ 


L+  -  L_ 


on  I+r 


with 
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*  ^  * 

and  L  £  H  . 
cause  of  the 
the  boundary 


=  Li  on  IMI  +  r)  u  !I+n)\(l  +  r) 

This  guarantees  that  |j  L  5 1  ^  "  L .  !|  while,  be¬ 
lt1  { I )  1 

exponential  decay  of  L  ,  "VL^  away  from 

of  I+r  . 


3.  CUBIC  SPLINE  INTERPOLATION  AT  KNOTS 

In  this  section,  we  establish  the  checkerboard  pattern 
for  A  in  case 

(7)  A  =  (Nj,4(tin 

of  cubic  spline  interpolation  at  knots. 

Given  that  the  interpolation  problem  is  correct,  we  now 
know  that  there  are  just  three  possibilities:  m+  =  0,1,2. 

Case  m+=0.  In  this  case,  VI  =  YL  ,  and  from  (4)  we  see  that 
A  1  is  upper  triangular.  Therefore,  as  discussed  in  Section  2, 
the  first  diagonal  of  A  is  main,  i.e.,  diagonal  1  in  our  way 
(7)  of  writing  A  .  It  follows  that 

(-)1+j+1A_1(i,j)  10,  all  i, j 

and,  in  particular,  the  solution  a  of  Aa  =  ( ( — ) 1 )  satis¬ 
fies  (-)i+1ai  >  0  ,  all  i.  Thus,  a  alternates  but  in  the 
wrong  sense  if  we  are  after  bounding  the  local  mesh  ratio  in 
terms  of  jj  R  ||  using  the  argument  of  Proposition  1.  In  fact, 
it  is  not  difficult  to  construct  a  knot  sequence  t  for  which 
the  interpolation  problem  is  correct  and  for  which  t^ — >0  as 
i  — -  -<*>  so  strongly  that  m  =  00  .  My  statement  at  the  last 
Bonn  conference  (  [4;p.48]  )  that  mfc  1  const  ||  R  jj  must  there¬ 
fore  be  qualified  to  exclude  the  case  m+  =  0  and  the  analogous 
case  m+  =  2. 

Case  m+=l.  We  find  it  convenient  to  associate  with  the  ele¬ 
ment  L  fell  the  nullspline 


L  :=  l  L ( j )  N. 

j  3 ' 

for  which  it  supplies  the  B-spline  coefficients.  We  claim  that 
(8)  L6tt'\{0!  implies  L*  (t^  L"  (t^  >  0  ,  all  i  . 

Indeed,  tl  contains  a  sequence  (L^j)  so  that 

L’li]  (ti>L(i]  (ti)  >  0  '  hence  L'[i]  {tj)L[i]  <tj)  >  °  f°r  3  -  if 


* 


i 


i 


I 
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I 


« 


I 


t 

I 


by  [1].  Since  Yl  is  finite-dimensional,  a  properly  normalized 
subsequence  then  has  a  limit  L  for  which  L’ (t^)L"  (t^)  0  , 

all  i  .  But  then,  by  (1], 

L*  (ti)L"  (t  )  >  0  ,  all  i 

!|L!|i+j  >  2jiiL|!.  ,  all  j  >  i 

with 

IN!  j  :=  max  {|L'  (t.)  !,  |L"(t.)  |}  . 

Suppose  now  that,  by  way  of  contradiction,  M  £  'll  '*{0}  satis¬ 
fies  M’  ( t ^ ) M " ( t . )  5  0  for  some  i  .  Then,  again  by  [1], 
M'(tj)M"(t.)  <  0  for  j  <  i  and 

II MH  i_ j  1  2j!iMl.  ,  j-1,2,3 . 

showing  that  then  M'  (t^)  and  M"(t^.)  both  would  increase 

exponentially  as  j - *  -*>  ,  while  M  itself  stays  bounded 

since  M«H  .  This  would  imply  that  At^  decreases  exponent¬ 
ially  as  j — «•  -*>  and  then,  given  that  2~3||L|li  >_  |[ L |(  ^ _ j  for 
j  =  l ,  2  , . . .  ,  L  would  surely  also  be  in  contradicting  the 
fact  that  m-  =  1  . 

One  proves  analogously  that 

(9)  LfeV-+s{0}  implies  L'  (t^)  L"  (t^)  <  0  ,  all  i  . 

(8)  and  (9)  imply  (see  [1])  that 

(10)  for  all  LfeflNfO),  L"  (t^  L"  (t±+1 )  <  0  ,  all  i  . 

Next,  we  claim  that 

(11)  for  all  L6H\(0),  L.  L"  ( t^+2)  <  0  ,  all  j  . 

For  this,  recall  (e.g.,  from[2;p.270] )  that  f  =  la.N.  .  im- 

3  3  ’ 4 

plies 

(12)  a  =  I  (-)3-r*(3-r)  (T)f (r)  (T) 

J  r<4 

for  any  xe(tj,t^+4),  with 

i(i(x)  :=  (tj+1-x)  (tj+2-x)  (tj  +  3-x)/3!  , 

Since  L  =  £  L .  N .  •  ,  we  then  get 

“J  J  •  q 

L.  =  *-(tJ+1)L'(tj+1)  -  ♦Mtj+1)L-(tj+1) 

while  ij/"(ts+1)  >  0  >  (tj  +  1)  .  Therefore,  with  (8), 
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0 


V(tj+ i>  = 


*"<V1)L,(V1,L"(V1)  ' 

^(tj+1)fL“(tj+1)]: 


if  L  £  ft  MO},  and  (10)  now  finishes  the  proof  of  (11)  for 
this  case.  The  case  L  €  ft  MO]  uses  (12)  with  x  =  t..+3 
stead . 


m- 


Finally,  let  now  L 
-1 , 


be  one  of  the  columns  of 


(*,i) 
+ 


L  =  A 

say.  Since  m_  =  m 

f  L“(j)  , 

(13)  L  ( j )  =  { 

l  L+  ( j  )  , 


1 ,  we  have 
for  j  <  i+2 
for  j  >  i 


with  neither  L  nor  L  just  0  .  Consequently,  by  (10)  and 
(11),  L  alternates  (strictly)  and,  as  to  the  sense  of  that 
alternation,  we  have  from  (12)  (with  j  =  i-2)  that 

A_1  (i-2,i)  =  <|»"  (ti_1)L'  (ti_1)  -  *•  (ti_1)L"(ti_1) 

while,  from  (13),  (3)  and  (9), 

4<"(ti_1)  >  0  >  'i>'(t._1)  ,  L*  (ti„1)L"  (ti_1)  >  0 

r(fci+1),  *'(ti+1)  <  0  ,  L'<ti+1)L»<t.+1)  <  o 

We  conclude  that 


sign  A  1(i-2,i)  =  sign  L" (ti_1)  =  sign  L"(ti+1)  . 

But  this  implies  that  A-^(i-2,i)  >  0  ,  since  the  contrary  as¬ 
sumption  would  give  L'  (tT)  >  0  >  L’  (tt)  (since  =  1) ,  an 

impossibility. 

In  conclusion,  if  m”  =  m  =  1  ,  then 
(14)  (-)i+^A  1 < i , j )  >  0  ,  all  i,j. 


and,  in  particular,  the  solution  of  the  linear  system  Aa  = 
((-)1)  does  satisfy  the  condition  (-)xa^  >  0  ,  all  i,  needed 
in  the  argument  for  the  finiteness  of  the  local  mesh  ratio. 

Note  that  we  proved  (14)  here  without  recourse  to  finite 
sections.  Even  in  this  simple  case,  I  still  do  not  know  whether 
the  matrix  A  has  a  main  diagonal  (though  I  don't  think  it 
would  be  very  hard  to  proved 
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